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O . We study multi-user transmission and detection schemes for a multi-access relay network (MARN) with 



o 
o 



linear constraints at all nodes. In a (J,J a ,R a ,M) MARN, J sources, each equipped with J a antennas, 
communicate to one A/-antenna destination through one i? a -antenna relay. A new protocol called IC-Relay- 



> 

TDMA is proposed which takes two phases. During the first phase, symbols of different sources are transmitted 

OO ■ 

concurrently to the relay. At the relay, interference cancellation (IC) techniques, previously proposed for systems 



with direct transmission, are applied to decouple the information of different sources without decoding. During 
the second phase, symbols of different sources are forwarded to the destination in a time division multi-access 
(TDMA) fashion. At the destination, the maximum-likelihood (ML) decoding is performed source-by-source. 

•rH . 

The protocol of IC-Relay-TDMA requires the number of relay antennas no less than the number of sources, 
i.e., R a > J. Through outage analysis, the achievable diversity gain of the proposed scheme is shown to 
be min{ J a (R a — J + l),Ra,M}- When M < J a (l — TT^)' ^ e P ro P osec l scheme achieves the maximum 
interference-free (int-free) diversity gain R a M. Since concurrent transmission is allowed during the first phase, 
compared to full TDMA transmission, the proposed scheme achieves the same diversity, but with a higher 
symbol rate. 

Index Terms: Multi-access relay network, distributed space-time coding, interference cancellation, 
orthogonal designs, quasi-orthogonal designs, cooperative diversity. 

Part of this work was presented at IEEE Wireless Communication & Networking Conference (WCNC) 2010. 
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I. Introduction 

Node cooperation improves the reliability and the capacity of wireless networks. Recently, many 
cooperative schemes have been proposed [l]-[4], and their multiplexing and diversity gains are 
analyzed. Most of the pioneer works on cooperative networks focus on cooperative relay designs 
without multi-user interference by assuming that there is one single transmission task or orthogonal 
channels are assigned to different transmission tasks, e.g., [l]-[4]. As a general network has multiple 
nodes each of which can be a data source or destination, multi-user transmission is a prominent 
problem in network communications. 

One model on multi-user cooperative communication is interference relay networks [5]. Multiple 
pairs of parallel communication flows are supported by a common set of relays. Each source targets at 
one distinct destination. Two transmission schemes using relays to resolve interference were proposed. 
The zero-forcing (ZF) relaying scheme uses scalar gain factors at relays to null out interference 
at undesired destinations [6]-[8]. The minimum mean square error (MMSE) relaying scheme also 
uses scalar gain factors at relays but to minimize the power of interference-plus-noise at undesired 
destinations [9], [10]. Both relaying schemes require the gain factors calculated at one centralized 
node having perfect and globe channel information, then fed back to the relays. These papers discuss 
the multiplexing gain and designs of the optimal scalar gain factors, but do not provide diversity gain 
analysis. In addition, for general multi-user cooperative networks, where communication flows may be 
unparallel, these schemes cannot be applied straightforwardly. For example, for a network in which 
several sources have independent information for one single-antenna destination, the ZF and MMSE 
relaying cannot resolve information collision at the destination. 

In this paper, we consider a multi-access relay network (MARN), in which J sources, each equipped 
with J a antennas, send independent information to one M-antenna destination through one i? a -antenna 
relay. We denote this network as a (J, J a ,R a ,M) MARN. For MARNs, a straightforward scheme is 
to use full time division multi-access (TDMA), where sources are allocated orthogonal channels for 
both hops of transmissions. Distributed space-time code (DSTC) [2], [11] is performed at the relay to 
gain high diversity performance without any channel state information (CSI). Such a scheme with full 
TDMA and DSTC at the relay is denoted as full-TDMA-DSTC. It achieves the maximum diversity 
gain R a mm{J a , M}. Since interference is avoided, this diversity gain is denoted as interference-free 
(int-free) diversity and provides a natural upperbound on the spatial diversity gain for all multi- 
user transmission schemes in the MARN. However, the spectrum efficiency of full-TDMA-DSTC is 
low. Another intuitive scheme is to allow multi-user concurrent transmission in both hops. The relay 
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conducts decode-and-forward (DF) by jointly recovering all sources' symbols. However, the decodings 
at the relay and the destination induce high processing complexity which is exponential in the number 
of sources. In [12], ZF beamformers are used in networks with two sources to make sources' signals 
orthogonal at the destination. The relay uses amplified-and-forward (AF). Nevertheless, the beamformer 
coefficients or global channel information need to be fed back to the sources, which induces a high 
protocol cost. In [13], we proposed a scheme called DSTC-ICRec that does not require CSI feedback 
to the relay and sources. The scheme allows concurrent transmission in both hops and uses multiple 
destination antennas to perform interference cancellation (IC) [14]— [16]. However, it trades overall 
diversity for spectral efficiency and cannot achieve the int-free diversity [17]. 

Based on the above discussion, a new scheme, called IC-Relay-TDMA, is proposed in this paper 
to allow multi-user concurrent transmission in the source-relay link. The multi-user interference is 
canceled at the multi-antenna relay by the linear IC techniques proposed in [14]— [16]. Then, space- 
time block code (STBC) and TDMA are used for the relay to forward signals of different sources to 
the destination. The merits of this scheme is summarized as follows: 

1) The IC-Relay-TDMA scheme applies not only to MARNs but also to general multi-user cooper- 
ative networks with multiple destinations and arbitrary patterns of communication flows as long 
as R a > J. The scheme requires CSI to be available at the receiving nodes only and no feedback 
is needed. The relay processing is linear and the decoding complexity at the destination is linear 
in the number of sources. 

2) It is proved rigorously that IC-Relay-TDMA achieves a diversity of min{ J a (R a — J + 1), R a M} 
in a (J, J a , R a , M) MARN. 

3) When M < J a ^1 — IC-Relay-TDMA achieves the int-free diversity, which is the max- 
imum spatial diversity achievable for (J, J a , R a , M) MARNs. Thus, the concurrent first-step 
transmission of the scheme induces no diversity penalty while improves the spectrum efficiency. 
The symbol rate of IC-Relay-TDMA is symbols/user/channel use where R denotes the 
symbol rate of the STBC used in the relay-destination link. Since the symbol rate of full- 
TDM A-DSTC is jjjTRo)' IC-Relay-TDMA achieves the same diversity, but with higher symbol 
rate, compared to full-TDMA-DSTC. 

The rest of the paper is organized as follows. Section II provides the network model. Section III 
introduces the IC-Relay-TDMA scheme. Its achievable diversity and symbol rate are discussed in 
Section IV. Section V shows numerical results and conclusions are given in Section VI. Involved 
proofs are presented in appendices. 
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Notation: For a matrix A, let A', A*, and A be the transpose, Hermitian, and conjugate of A, 
respectively. || A|| is the Frobenius norm of A. tr {A} calculates the trace of A. ® denotes Kronecker 
product. I n is the n x n identity matrix. 0„ is the n x n matrix of all zeros. For two matrices of the 
same dimension, A y B means that A — B is positive definite. f(x) = o(x) means lim + ^^ = 0. 
E [x] denotes the expected value of the random variable x. 

II. Network Models 

In this section, we describe two network models that are used in the paper. First, we introduce the 
MARN, then the general multi-user cooperative network. 

Consider a relay network with J sources each with J a antennas, one relay with R a antennas, and one 
destination with M antennas. There is no direct connection from sources to the destination, because 
the sources are far from the destination. The system diagram is shown in Fig. 1. 

Denote the fading coefficient from Antenna k (k = 1, . . . , J„) of Source j (j — 1, . . . , J) to Antenna 
% (i — 1, . . . , R a ) of the relay as fj^\ The J„xl channel vector from Source j to relay Antenna % is 
denoted as whose A;-th entry is fjf). Denote the fading coefficient from relay Antenna i to Antenna 
m of the destination (m = 1, . . . , M) as gi m . The R a x 1 channel vector from the relay to destination 
Antenna m is denoted as g m whose i-th entry is g im . All fading coefficients are assumed to be 
identically and independently distributed (i.i.d.) with CA/"(0, 1) distribution. We assume a block-fading 
model with coherence interval T. 

To allow IC at the relay, we assume that R a > J. This can be realized through user admission control 
in the upper layers. We assume that sources have equal numbers of antennas. Our proposed protocol can 
be extended straightforwardly to networks where sources have unequal number of antennas. Further, 
to focus on the diversity performance, all sources and the relay are assumed to have the same average 
power constraint P. The extension to nonuniform power constraint is also straightforward. Throughout 
the paper, we assume global CSI at the destination; but the relay has only the backward CSI, i.e., 
channel information from sources to the relay. The channel information can be obtained by sending 
pilot sequences from sources and the relay [11], [18]. No feedback or channel estimation forwarding 
is required. Perfect synchronization at the symbol level is assumed for the network. 

For complexity considerations, two linear constraints are imposed on the network. For one, the relay 
linearly transforms its received signal to generate its output signal without decoding. For the other, 
the decoding complexity at the destination is linear in the number of sources. It can be verified that 
the full-TDMA-DSTC scheme mentioned in the introduction section and the DSTC-ICRec scheme 
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proposed in [13] satisfy these two linear constraints. 

The general multi-user cooperative network has J + N + 1 nodes. J multi-antenna sources, denoted 
as S = (si, s 2 , . . • , s,j), send independent information to iV multi-antenna destinations, denoted as 
D = (di,rf 2 , • • • , d N ), through one i? a -antenna relay. In an indoor environment, the mobile stations 
can be the sources and the destinations, and the access points connecting through cables can be the 
relay. The set of sources from which Destination n receives information is described as I n = {s, s e S}. 
The profile of communication flows of the whole network can be described as / = I 2 , ■ ■ ■ , In) - For 
example, a network with three single- antenna sources S = (si, s 2 , s 3 ), two single- antenna destinations 
D = (di,d 2 ), and one four-antenna relay is shown in Fig. 2. Destination 1 receives information 
from Sources 1 and 2, while Destination 2 receives information from Sources 2 and 3. The profile 
of communication flows can be expressed as / = ({s 1: s 2 }, {s 2 , s 3 }). Specifically, when J = N 
and I = ({«i}, {s 2 }, • • • ,{sat}), the network becomes the interference relay network with parallel 
communication flows. When N = 1 and / = ({si, • • • , sj}), the network becomes the MARN. All 
nodes are assumed to be synchronized. Extension to asynchronous networks is straightforward using 
the random access and IC methods in [19]. Although for the clarity of presentation, we present our 
protocol using the MARN model, we will show that it can be applied straightforwardly to this general 
network model. 

III. IC at the Relay: IC-Relay-TDMA 

It is well known that for cooperative networks relaying can improve communication reliability and 
coverage. In this paper, we show a new dimension in the design of multi-user relay networks: IC at 
relays. In our MARNs, to improve the spectral efficiency, we allow concurrent multi-user transmission 
in the link between the source and the relay. Since this source-relay link is a multi-antenna multi- 
access channel (MAC), the multi-antenna relay has the potential to cancel the induced multi-user 
interference. Cancelling interference at the relay improves the signal to interference-plus-noise ratio 
of the relay-destination link and simplify the signal processing at the destination. Thus, this idea has 
the potential of providing good performance when the relay-destination link is the bottleneck of the 
network. Based on the above considerations, we propose a protocol called IC-Relay-TDMA, in which 
the relay conducts IC using linear transformations but not decoding before forwarding int-free signal 
to the destination by TDMA. In Subsection III-A, we describe the protocol for general (J, J a ,R a , M) 
MARNs. Its application in a general multi-user cooperative network is discussed in Subsection III- 
B. Then, the use of IC-Relay-TDMA in one simple network is illustrated in Subsection III-C as an 
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example. 



A. The Protocol of IC-Relay-TDMA 

In this subsection, we explain the protocol of IC-Relay-TDMA. The protocol consists of two phases. 
During the first phase, all sources send information to the relay simultaneously using STBC with 
ABBA structure [20], [21]. The relay overhears superimposed signals of all source information and 
conducts multi-user IC [16]. During the second phase, the relay conducts a scheme called MRC- 
STBC to enhance the communication reliability in the second transmission, then forwards information 
of different sources in TDMA to the destination. The destination decodes each source's information 
independently. A block diagram of IC-Relay-TDMA is shown in Fig. 3. The details of the protocol 
and corresponding formulas are described in the following. First, we consider the scenarios when J a 
is a power-of-2, then extend to the cases that J a is not a power-of-2. 

Phase 1: When the number of antennas at each source is a power-of-2, i.e., J a = 2 n , n e J\f, each 
transmitter constructs 4f constellations (e.g. PSK, QAM constellation and their rotations), denoted as 
S u , u — 1, . . . , 4f. The average power of the constellations is normalized to be one. The constellations 
need to satisfy the following condition for diversity gain: 



One approach to construct such constellations is through rotation [22], [23]. For example, when J a = 4, 
i.e., n = 2, two BPSK constellations can be constructed as Si = { — 1,1} and S 2 = {j, — j}, where S 2 
is rotated from Si by n/2. It can be verified that Si and S 2 satisfy the condition in (1). 

Source j independently and uniformly collects J a symbols s^\k = l,...,J a from these J a /2 
constellations with s2-i>4« e «S U . Then, a J a x J a STBC with ABBA structure [20], [21] is 
constructed by 



where the function S n maps J a variables to a J a x J a matrix through an iterative procedure as 



with Si(si, s 2 ) an 2 x 2 Alamouti code. All sources transmit simultaneously in this phase. The length 
of this phase is Ti = J a time slots. It is thus assumed that the coherence interval T is no less than J a . 




(1) 
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Relay Antenna i overhears a 7\ x 1 signal vector as 

3=1 * 

where Vj denotes the J a x 1 additive white Gaussian noise (AWGN) vector, whose r-th entry v Ti is 
i.i.d. CJ\f(0, 1) distributed. Note that the first phase transmission is virtually a multi-antenna multi- 
access channel. When R a > J, the IC scheme originally proposed for direction transmission [16] can 
be conducted at the relay to fully cancel the multi-user interference. In [16], the IC procedure was 
discussed explicitly only for a system with at most four antennas at each source. Here, we describe 
this procedure for a general system with J a = 2 n antennas at each source. Without loss of generality, 
we discuss how the relay cancels interference from Source 2 to Source J and obtains an int-free 
observation of the information of Source 1. 

The IC procedure has two steps. In the first step, the relay separates the system that communicates 
J a = 2 n symbols for each source into 2 n ~ 1 equivalent Alamouti systems. In the second step, the 
IC scheme in [15] is applied to each Alamouti system to iteratively cancel interference from Source 



2 to Source J. Denote h/ as the l-th row of the 2" 1 x 2" 1 Hadamard matrix H 



*1 7*3 ' 
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„Ci) Jj) 
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U) 



1 1 



f 

i to 



J li ' J 3i ' " " " ' J I 



it 



(2»-l)i 



f 



(3) _ 



n -i. Let s, 
1 1 



U) 



J2i > Jii i ■ ■ • J J2 n i 



[vii,v 3i , . . . ,W( 2 n-i)i]*, and v ic = [v 2 i, v Al , . . . , U( 2 »)j] As the first step, relay Antenna % calculates r H 
(hi <g) I 2 )rj to obtain equivalent Alamouti systems as follows, 



h)Ti 



j=l:J 



+ 



, 1 = 1,.. .,2 



n-l 



(3) 



Denote the first and second entries of ru as fm and r^ 2 , respectively. Due to the Alamouti structure 
of Si (hiSo \ h/Sc^j. Eq. (3) can be equivalently rewritten as 












+ 




-h;v ie 



(4) 



3 C#) 



— 



For the second step, the relay cancels interference for each Alamouti system in multiple iterations, 
where interference of one source is cancelled in each iteration. Stack r H and at different relay 



antenna as f; 



n Ra ] 1 and F} 



(.3) 



? (3)t 
11 ' 



F (3)t 



. Denote Ti(i) as the 2(R a 



l)x 



2(R a — i) IC matrix to cancel Source J — i for System /; rj(i) and Fi(i) as the remaining 2(R a — i) x 1 
signal vector and the remaining 2(R a — i) x 2J equivalent channel matrix after cancelling Source 
J — i + 1 for System respectively. The iterative IC procedures are as follows: 
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. Initialization: F,(0) = [fJ 1} , . . . , ] 
• Iteration: For % = to J — 2, 

1) Form the IC matrix as 



rj(0) = r,. 



2F, (0 ; - , 

l F (i) ||2 

2F, (!) ; - . 

I, J— i,l 
;|F« ||2 



\\*%-i,2 II 2 

o 2 



2F 



o 2 

(i). 



T (0 II 
■ J,J-i,3ll 



2F 



(i). 

i,J-i,l 



o. 



o 2 



2F 



o 2 
o 2 

(i). 

l,J-i,M- 



■ l,J-i,M~ 



(5) 



-I 2(.R -i-l)x2(.R -i) 



2) Cancel interference of Source J — iby multiplying r ; (i) with J-^i). The remained equivalent 
received signal can be calculated as r ; (i+l) = Fi(i)ri(i) and the remained equivalent channel 
matrix can be calculated as F t (i + 1) = jFi(i)Fi(i). 
The 2x2 matrix F^ in (5) denotes the (p, g)th 2x2 submatrix of F t (i). After J - 1 iterations, 
information of Sources J to 2 is cancelled and the remaining signal vector r;(J — 1) only contains 
information of Source 1. The overall IC matrix that jointly cancels Sources 2 to J is T\ = Yli=o -^KO- 
To help the presentation, let = r ; (J — 1). From this iterative procedure, we have 



f, (1) =^« = J^«F I (1) S l (1) +7iv,, 

V •-'a 



(6) 



where v ; = [v* 1? v* 2 , 



A 2(R a — J + 1) x 1 vector observation of Source l's information 



is carried in r|^. Eq. (6) implies that the rows of T\ are in the null spaces of the columns of F^ to 
F\ J \ i.e., J-iF^ = for j = 2, . . . , J. Thus, this IC process is an iterative realization of ZF. Different 
from conventional ZF which uses pseudo-inverse of the channel matrix, the IC method needs no 
channel matrix inversion. Similarly, the relay can obtain int-free vector observations of other sources' 
information. 

Phase 2: In this phase, the relay conducts a process called MRC-STBC [24], then forwards 
information of each source to the destination in different time slots. The destination decodes source- 

(i) 

by-source and jointly recovers the symbols contained in So . Without loss of generality, we only 
consider how Source l's information is processed by the relay and decoded at the destination. The 
maximum-ratio combining (MRC) step is first conducted to maximize the SNR of rj 1} . The MRC can 
be represented as 



s(i) 



2F( 1 )* < 7>( < 7y>)-i ? ! 



0-h 




^ + 



2F ( [ 1) *f;(f 1 f;)- 1 f 1 v 1 



0-h 



(7) 



-,(!) 
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The entries of the vector after MRC, are soft estimates of the entries of s|^. From (3), (4), (6), and 



- (1 u — : Ja : : mr^2, which implies that 



(7), the covariance matrix of can be calculated as 

the two noise elements in vj 1 ^ are i.i.d.. Also, the noise vectors vj 1 ^ of different Alamouti systems are 
independent, due to the orthogonality of the Hadamard matrix H n _i, but not identical. Following the 
MRC step, to forward Source l's information, the relay uses generalized orthogonal STBCs to encode 
entries of . We especially consider generalized orthogonal STBCs because of its full diversity and 
symbol-wise decoding [25]. Other designs such as quasi-orthogonal STBCs [21] can also be applied but 
with higher decoding complexity. In general, consider using a T 2 x R a generalized complex orthogonal 
design that carries K symbols. Note that each contains information of two symbols. The relay 
waits for K symbols from \K/2] Alamouti systems, denoted as sf \ k = 1, . . . , K (the subscript I is 
removed without confusion for conciseness), to generate the T 2 x R a codeword as 



(8) 



k=l:K 



where tj is the T 2 x 1 signal vector to be transmitted at relay Antenna i; A k and B fc are T 2 x 
R a relay encoding matrices for generalized orthogonal designs as found in (4.67) in [26]; and c = 
. / D ^ P PT2 i , is the power normalization coefficient at the relay. Since the processing in (3), (6), 

V a y~ + 2R a -2.I + l I 

(7), and (8) are linear, the transmitted signal vectors at the relay are linear in its received signal vectors 
i-j. Assume that the coherence interval T is no less than T 2 . The relay concurrently forwards tj on 
Antenna i. The received T 2 x 1 signal vector at destination Antenna m can be expressed as 

E (Re{4 1) }A fc +jIm{4 1) }B A; )g m + u m , 



X m = E **im 

i=l:R a ' k=l:K 

where w m denotes the T 2 x 1 AWGN vector at destination Antenna m; u m denotes the equivalent noise 

vector, u m = c (R- e {^^l-^-fc + j im {^jBfc ) g m + w m with the additive noise defined in 

k=l:K \ " ' 



+ w m = c 




(9) 



(7). 

At the destination, a 2T 2 M x 1 vector is formed by stacking x m into x = [Re {xi}*, Im {xi}*, 
Re {x M }', Im {xm} 1 ]'. After straightforward calculation, the system equation can be written as 
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where 





Re{g m } 










Im{g m } 





-Im{B fe }-Re{B fc } 
Re{B fc } -Im{B fc } 



Re{A fc } -Im{A fc } 
Im{A fc } Re{A fc } 

With generalized complex orthogonal designs, G*G = ||G|| 2 I 2 a-, where G = [gi, • • • , gjw]. Denote 
a fc as a 1 x 2K vector whose 2k — 1 and 2k entries are 1 and j, respectively, and all the other entries 
are zero. The destination can obtain a soft estimate of by the following calculation, 

G x / P 



IGIP ~ C V J n Sk +CVk ' IIG 



(ID 



where w k is the equivalent noise with CAf(0, 1) distribution and independent for different k. x k is 
a soft estimate to which is a linear superposition of Source l's symbols from (4). Without 



loss of generality, we assume that provides a soft estimate to h. k s ( L> ,k = 1,...,2™ 1 . Denote 



(i) 



Xi ... X 2 n-i 



v(l) = 



,-,(!) .-.(1) 



. V 



2 n-l 



and w = 



Wl . . . 1W2™- 1 



From (11), we have 



x 




H n _ lS « + cvw + 



(i) 



w 



G 



(12) 



The destination performs ML decoding to decode based on (12) as 



arg mm 

s 




S.T 1 



x — c 




H„_is 



(13) 



where S u is the covariance matrix of the equivalent noise vector u. After straightforward calculation, 
we have 



W-l„2 



c diag 



1 1 

7i ' 72 ' 



72™ -i 



H 12"- 1 ; 

7 9 



(14) 



where 7, = |tr {f{ 1} *^* (J 7 ^* )- 1 ^fJ 1) } and lg = ||G|| 2 . Similarly, s^ 1} can be jointly decoded. 
Transmission of other sources' symbols can be performed similarly using orthogonal time slots. To 
decode all symbols from all J sources, the destination only needs to conduct 2 J procedures of ML 
decoding of 2 n_1 symbols. The complexity is linear in the number of sources. 

When J a is not a power-of-2, 2™ x 2 n quasi-orthogonal STBCs with ABBA structure are used with 
2 n the smallest power-of-2 number greater than J a . During the first phase, each source concurrently 
transmits the first J a columns of the block codes in Ti = 2 n time slots. Similar to the case when J a is 
a power-of-2, the resulting multi-user interference can be cancelled using Eq. (3), (4), and (6) at the 
relay by treating fjfj = for k = J a + 1, . . . , 2 n . During the second step, symbols of different sources 
are forwarded by MRC-STBC in TDMA. Symbols are decoded source-by- source at the destination. 
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B. Application in General Multi-User Cooperative Networks 

IC-Relay-TDMA can be applied to the general multi-user cooperative networks with multiple 
destinations. During the first phase, all J sources send information to the relay concurrently. The 
relay separates multi-user signals using IC without decoding. During the second phase, int-free soft 
estimates of each source's symbols are encoded using STBC. The relay broadcasts each source's block 
codes using TDMA. All destinations receive int-free signals from all sources. The destination decodes 
its desired information and discards undesired information. Note that transmission and processing at 
the relay do not depend on the relay-destination link and the number of destinations. IC-Relay-TDMA 
is robust to the dynamics of the destinations and no relay-destination channel information is required 
at the relay. On the contrary, for ZF relaying and MMSE relaying, the relay needs to acquire the 
channel information of new destinations and updates the scalar gain factors, which takes substantial 
protocol overhead. IC-Relay-TDMA can be applied to any patterns of communication flows when 
R a > J, but ZF relaying and MMSE relaying require the flows to be parallel. It should be mentioned 
that IC-Relay-TDMA has a lower symbol rate than that of ZF and MMSE relaying. For the same bit 
rate, larger constellations are required. 



C. An Example: IC-Relay-TDMA for a (2, 2, 2, 1) MARN 

In this subsection, we present one example of using IC-Relay-TDMA in a (2, 2, 2, 1) MARN, where 
there are two double-antenna sources, one double-antenna relay, and one single- antenna destination. 
The description of the proposed scheme in the previous subsection is lengthy as it is for a general 
MARN setting. In this network example, we will see that some processing are naturally simplified or 
become unnecessary and the main ideas behind the scheme are more clearly illustrated. The complexity 
at the relay and the destination can be further reduced. 

During the first phase, only one constellation is required and the constraint on the constellation in (1) 
becomes trivial. Both sources collect two symbols from the same constellation, and concurrently send 
two Alamouti codes, i.e., S^ = Si (s^\s{^ , j = 1,2. Since there is one Alamouti system only, 
the signal separation illustrated in (3) is also not needed. At the relay, only one round of IC iteration 

„p(2)« 2 p( 2 )* 

is needed. The interference of Source 2 can be cancelled by using the IC matrix T = <h— , — m — , 

j & y IIF< 2) p' ||F< 2 »pJ' 

where F-^ is the Alamouti channel matrix from Source j to relay Antenna i, i.e., F-^ = Si (^fij\ f^ij ■ 
A 2 x 1 vector observation of Source l's symbols can be obtained from (6). During the second phase, 
after the MRC represented in (7), the vector containing information of each source is encoded into an 
Alamouti block code and forwarded to the destination in TDMA. At the destination, the equivalent 
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system equation for Source 1 can be written as 



Xi 



<-\l g-Si (gi,g 2 ) 



S (!) 



(15) 



u 



where x T , v^, and w denote the received signal at time slot r, the equivalent noise vector of the relay 
in (7), and the 2 x 1 equivalent noise vector at the destination, respectively. The covariance matrices 
of v* 1 ) and w are 7/I 2 and I 2 , respectively, where 7/ = itrjF^^J 7 * (FF*)^ 1 FF^}. Again, for 
this simple network, the processing in (10) is not needed. The destination directly performs the ML 
decoding based on (15), which, for this network, simplifies to 



argmm | x - c\j ^Si (gi, g 2 ) s (1) ) S u : ( x - c\j ^Si (gi, g 2 ) s (1) ) , 



(16) 



2 J \ V 2 

where S u is the covariance matrix of the equivalent noise u and after straightforward calculation, 
E u = (c 2 7/ (|(?i| 2 + Ifi^l 2 ) + l) 12- Since S u is a multiple of identity independent of the information 
symbol, it can be omitted in the ML decoding without any performance loss. Due to the Alamouti 
structure in Si (gi,g 2 ), (16) can be further decomposed into two symbol-wise decodings. 



IV. Performance Analysis 

In this section, we provide diversity analysis for the protocol of IC-Relay-TDMA and discuss its 
properties. Subsection IV-A is on the diversity analysis. In Subsection IV-B, we discuss the symbol 
rate of the scheme and when it achieves the int-free diversity. 



A. Diversity Analysis 

The diversity of a communication system is defined as the negative of the asymptotical slope of the 
bit error rate (BER), d = — Jim l ° o | p . For fixed constellations, can be replaced with pairwise symbol 
error rate (SER). Since the ML decoding in (13) is source-by-source and the network parameters and 
processing at the relay and the destination are statistically homogenous, the diversity of each source 
is identical. We only need to analyze the diversity of one source, without loss of generality, Source 
1. The concatenation of two hops of transmission and relay processing make the calculation of SER 
extremely difficult. Thus, to aid the diversity gain analysis, in the following lemma, we provide a 
method to calculate the diversity based on a formula with the outage probability structure without 
explicitly calculating the SER. 
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Lemma 1: Define the instantaneous normalized receive SNR as 7 = tr {H*_ 1 S~ 1 H„_ 1 }. When 

the constellations satisfy (1), the diversity of the ML decoding in (13) can be calculated as 

rf=1 . m logP( 7 < f ) 
e-^o+ log e 

Proof: See the appendix for the proof. ■ 
Lemma 1 says that diversity can be calculated using the outage probability of the instantaneous 
normalized receive SNR. Thus, diversity can be obtained from the minimum exponent of P(7 < e). 
More precisely, a random variable 7 provides diversity d if P(7 < e) = ce d + o(e d ) where c is a 
constant independent of e. Before the diversity theorem, the following lemma is introduced. 

Lemma 2: Let 71, 72, ... , 7tv, 7 9 be N+l instantaneous normalized receive SNRs. 7 3 is independent 
of 7 n for n — 1, 2, . . . , N. 7 9 provides diversity d\\ 7„ provides diversity d 2 . If 7 = 7 "7 9 , 

n=l:N ' n=l:JV 7 " 79 

7 provides diversity min-^, d 2 }. 

Proof: It can be shown by straightforward calculation that 7n7g < min{ 7„, N^}. The 

„=l:Ar 7n+79 n =l:N 

right-hand side has diversity min{rf 2 , d{\ from Lemma 1. Therefore, the diversity of 7 is upperbounded 
by min{<i 2 , d{\. To show the lowerbounds on the diversity, the following events are defined: 8 = {71 < 

72 < . . . < 7a?}, £ n = £ D{7n-i < 1 9 < In} for n = 1, . . . , N + 1 where 70 = and 7^+1 = 00. 
Since {£ 1 ,£ 2 , . . . , £n+i} is a partition of E, we have 

P( 7 <e,£) = P(7<e,5i) + P(7<e,5 2 ) + ... + P( 7 <e,5 JV+ i) 

< (P( 7 < e|5i) + P( 7 < e|£ 2 ) + . . . + P( 7 < e|W)) max{P(£ n )} 

< (^P( 7fl < 2e) + P ^ 7 n < 2ejj max{P(£ n )} = ce d + o(e d ), 

where d = mm.{di,d 2 } and c is a constant independent of e. For the third inequality, we have used 
the facts that P (7 < e\£ n ) < P (^p^ < e\£^j < P{^ g < 2e) when n < N, and P (7 < e|£jv+i) < 



P I S 7« < 2e ) • In the third line, the term max{P(£ n )} is independent of e, hence does not 

\n=l:N J 

affect the diversity. This is true for any orders of the sequence 71, ... ,7^. Note that P(7 < e) = 
2~2P(l < e, £) where the summation is over all possible orders. The diversity of 7 is lowerbounded 

£ 

by the minimum of the exponents of P(7 < e, £), which is min{rf 2 , d{\. Therefore, the diversity is 
min{c?2, di}. ■ 
Theorem 1: In (J, J a , R a , M) MARNs, IC-Relay-TDMA achieves a diversity of min{J a (P a - J + 
l),P a M} when R a > J. 
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Proof: From (13) and (14), the instantaneous normalized receive SNR can be calculated as 

7 = tr {H*_ 1 S~ 1 H n _i} = £ + (18) 

Since 7; is identical to the instantaneous normalized receive SNR in a multi-antenna multi-user 

i=l:2™- 1 

system with J J a -antenna users and IC at the iV-antenna receiver, 2n \ c2 1i provides diversity 

1=1:2™— 1 

J a (R a — J + 1) [27]. Since 7 9 is a Gamma distributed random variable with degree R a M, 7 9 provides 
diversity R a M. By Lemma 2, 7 has diversity min{ J a (R a — J + 1), i? a ^T}- ■ 



B. Performance Discussions 

This subsection discusses the condition for the proposed scheme to achieve the int-free diversity, 
the symbol rate, and the complexity of the proposed scheme. Comparisons with other schemes are 
also provided. 

The int-free diversity condition: Theorem 1 says that IC-Relay-TDMA achieves diversity min{ J a (R a - 
J + \),R a M}. Recall that the int-free diversity is defined as the maximum achievable diversity for 
(J, J a , R a , M) MARNs without interference, which is R a min{ J a , M}. When 

M < J a (l - , (19) 

IC-Relay-TDMA achieves diversity R a M, which is equal to the int-free diversity under (19). Eq. (19) 
is then called the int-free diversity condition. This condition implies that M < J a , i.e., there are 
more independent paths in the source-relay link than the relay-destination link. For these networks, 
the bottleneck of transmission is the relay-destination link. Intuitively, when the source-relay link has 
extra degrees of freedom, they can be used for IC without degrading the total diversity. This is the 
basic idea behind IC-Relay-TDMA. Some examples of networks achieving the int-free diversity are 
(2, 2, 2, 1); (2, 4, 2, 1); and (2, 2, 4, 2) MARNs. To the best of our knowledge, in multi-antenna MAC, 
there is no IC method that achieves full single-user diversity. For MARNs, this is possible due to the 
extra relaying step. For networks satisfying (19), the source-relay link provides enough extra degrees 
of freedom to cancel interference at the relay. 

The symbol rate: During the first phase, each source sends Ti = 2 n symbols during Ti time slots. 
During the second phase, assume that the relay uses generalized orthogonal STBCs of dimension 
T 2 x R a to carry K symbols. Then, the total number of time slots in the second phase is JT ^ 2 ■ Let 
the symbol rate of the STBC code used in the second hop as R = ||. Thus, the symbol rate of each 
source is R — T x j (Ti + JT ^ 2 ) = j+° R ■ If rate-1 codes (e.g., Alamouti code) are used in the second 
transmission phase, we have R a — 1 and the symbol rate of the scheme is j^j. 

DRAFT 



15 

Complexity: We discuss the complexities in terms of the number of sources at the relay and the 
destination. Note that the relay needs to cancel the interfering signals from J — 1 sources to obtain 
the int-free signal from one source and there are J sources needed to be decoupled. The complexity 
of the IC is quadratic in the number of sources at the relay. At the destination, the complexity of ML 
decoding is linear in the number sources. Therefore, the relay has higher order of complexity than the 
destination in terms of the number of sources. 

Comparison with other schemes: We now compare IC-Relay-TDMA with other schemes in MARNs. 
Recall that the proposed IC-Relay-TDMA scheme has concurrent transmission in the source-relay link 
only. We first compare it with full-TDMA-DSTC, which uses TDMA to avoid multi-user interference 
in both hops. The second compared scheme is DSTC-ICRec [13], which allows multi-user concurrent 
transmission in both hops and IC at the destination to decouple signals of different sources. Finally, 
we introduce DSTC joint-user ML decoding, which is similar to DSTC-ICRec excepts that, instead of 
conducting IC then decoding each source's messages independently, the destination jointly decodes all 
sources' messages. Note that the decoding complexity of this scheme is exponential in the number of 
sources. Thus, DSTC joint-user ML decoding does not satisfy the linear constraint at the destination 
mentioned in Section II, but the other three schemes satisfy the linear constraints both at the relay and 
the destination. We compare diversity, symbol rates, and other properties of these schemes in Table I. 
The details on the diversity results in this table can be found in [28]. 

For networks satisfying the int-free diversity condition, IC-Relay-TDMA achieves the same diversity 
as full-TDMA-DSTC with a higher transmission rate. This is due to the concurrent transmission and 
diversity redundancy in the source-relay link. Though the symbol rate of DSTC-ICRec is higher than 
that of IC-Relay-TDMA, a higher dimension constellation can be used for IC-Relay-TDMA to achieve 
the same bit rate with faster decaying error probability. DSTC joint-user ML decoding achieves the 
maximum int-free diversity with a symbol rate higher than that of IC-Relay-TDMA. However, the 
decoding complexity of the DSTC joint-user ML decoding is exponential in the number of sources, 
which is very demanding when J is large. This implies that the proposed IC-Relay-TDMA trades 
decoding complexity for symbol rate without losing diversity for networks satisfying the int-free 
diversity condition. It should also be noted that IC-Relay-TDMA requires backward CSI at the relay, 
while the other three schemes do not require any CSI at the relay. Backward CSI can be obtained via 
training and does not need any feedback. 
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V. Numerical Results 

In this section, we show simulated BER performance of IC-Relay-TDMA and its comparison with 
other schemes. In all figures, the horizontal axis represents the average transmit SNR, measured in 
dB. Since the noises at all nodes are normalized, P is equal to the average transmit SNR. The vertical 
axis represents the BER. 

Fig. 4 is on the BERs of IC-Relay-TDMA under four network scenarios: Network 1: (2,1,2,1) 
MARN; Network 2: (2,2,2,1) MARN; Network 3: (2,4,2,1) MARN; and Network 4: (2,2,4,1) 
MARN. For the first three networks, Alamouti codes are used; and for Network 4, the rate 3/4 
generalized orthogonal STBC shown in (4.103) of [26] is used. All networks apply BPSK modulation. 
Networks 2, 3, and 4 satisfy the int-free diversity. Fig. 4 shows that for these three networks, the 
proposed scheme achieves the int-free diversity 2 for Networks 2 and 3; and 4 for Network 4. Network 
1 does not satisfy the int-free condition in (19). Fig. 4 shows that Network 1 has diversity 1, which 
is less than 2, the int-free diversity. These simulation results verify Theorem 1 . 

In what follows, we compare the proposed IC-Relay-TDMA (Scheme 1) with DSTC-ICRec (Scheme 
2), full-TDMA-DSTC (Scheme 3), full-TDMA-DSTC CIR (Scheme 4), DSTC joint-user ML decoding 
(Scheme 5), IC-Relay-TDMA DF (Scheme 6), joint-DF-TDMA (Scheme 7) in the (2, 1, 2, 2) and 
(2,2,2, 1) MARNs. Schemes 1, 2, 3, and 5 are discussed in Subsection IV-B. To rule out the effect 
of the difference in the CSI requirements for Scheme 1 (the relay needs to know its channels with the 
transmitters) and Scheme 3 (no channel information at the relay), Scheme 4, originally proposed in 
[29] for single-user relay networks, is included as well. In this scheme, the relay uses its knowledge 
of the backward CSI to equalize the phase shift of the source-relay link, then forwards information to 
the destination by Alamouti DSTC. To allow decoding at the relay, Scheme 6 is introduced, which is 
similar to Scheme 1 except that all sources' symbols are decoded after IC at the relay and re-modulated 
by the same constellation before forwarding. For Scheme 7, the relay jointly decodes symbols from 
both sources without IC before forwarding each source's information using Alamouti DSTC in TDMA. 
Note that Schemes 1, 2, 3, 4 satisfy the linear constraints introduced in Section II; whereas Schemes 
5,6,7 do not. For Scheme 5, the decoding complexity is exponential in J at the destination. The relay's 
decoding complexity for Schemes 6 and 7 are linear and exponential in J, respectively. Schemes 1, 
4, 6, 7 require backward CSI at the relay, but the other three schemes need no CSI at the relay. To 
achieve 1 bit/user/channel use for all schemes, the modulation constellations used for the schemes 
are 8PSK, QPSK, 16PSK, 16PSK, QPSK, 8PSK, 8PSK, respectively. Since the destination in the 
(2, 2, 2, 1) MARN has only single- antenna, Scheme 2 is excluded from the comparison in Fig. 6. 

DRAFT 



17 

We first compare the BER of Scheme 1 with the other linear schemes (Schemes 2, 3, 4). In the 
(2, 1, 2, 2) MARN (Fig. 5), Schemes 3 and 4 achieve a diversity gain of 2, which is the int-free diversity 
gain. Schemes 1 and 2 achieve a diversity gain of 1 only. For this network, since the int-free condition 
is not satisfied, the proposed Scheme 1 performs worse than Scheme 4 for the entire simulated SNR 
range. Although it is worse than Scheme 3 for SNR larger than 27 dB due to the diversity loss, it 
outperforms Scheme 3 when the SNR is smaller than 27 dB due to its higher symbol rate. Compared 
with Scheme 2, the proposed scheme is superior for the simulated SNR range. At the BER level of 
1(T 2 , it is about 10 dB better. For the (2, 2, 2, 1) MARN (Fig. 6), the int-free condition is satisfied. 
The figure shows that Schemes 1 and 4 achieve a diversity gain of 2, while the diversity of Scheme 3 
is slightly less than 2. This is because for Scheme 3, there is a logP factor in the error rate formula 1 . 
As P increases, the diversity should approach 2. For Schemes 1 and 4, the MRC and equalization 
eliminate the logP factor. Scheme 2 cannot be used for this network because the destination has 
only one antenna and cannot conduct full IC. The array gain of Scheme 1 is higher compared to 
both Schemes 3 and 4, since a lower-dimension constellation is used to achieve the same bit rate. 
At the BER level of 1CT 3 , it is better than Schemes 3 and 4 by 10 dB and 5 dB, respectively. From 
the comparison, we can conclude that Scheme 1 is expected to outperform other linear schemes for 
MARNs satisfying the int-free condition, e.g., the (2,2,2, 1) MARN. 

Next, we compare Scheme 1 with the schemes not satisfying the linear constraints (Schemes 5, 6, 
7). Scheme 1 is first compared with Scheme 5. Scheme 5 achieves the int-free diversity from Table I, 
thus naturally having better BER than Scheme 1 in the high SNR regime for the (2, 1,2,2) MARN 
(Fig. 5). For SNR smaller than 20 dB, Scheme 1 outperforms Scheme 5. In the (2, 2, 2, 1) MARN 
(Fig. 6), where both schemes achieve the int-free diversity, Scheme 1 outperforms Scheme 5 in the 
entire SNR, e.g., it is about 6 dB better when BER= 10 3 . The gain is obtained because in Scheme 1, 
user interference is avoided in the second step and the received signal quality is high. Then, Scheme 
1 is compared with Scheme 6, which has additional decoding after IC. From both Fig. 5 and 6, 
we can observe that there is no diversity improvement by additional decoding after IC for Scheme 
6. For the array gain, Scheme 6 outperforms in the low SNR regime (about 0.5 dB in Fig. 5 and 
1.3 dB in Fig. 6); and has the same BER as Scheme 1 in the high SNR regime. This is because the 
BER performance is mainly restricted by interference in the high SNR regime. Finally, Scheme 1 is 
compared with Scheme 7, which allows joint decoding at the relay. We can see that joint decoding at 

'if quasi-orthogonal designs are used as the distributed STBC, the logP factor does not appear and diversity 2 can be achieved as 
proved in [29]. However, the use of quasi-orthogonal designs requires the coherent interval to be 4. In this simulation, T = 2 and 
Alamouti designs are used at both the relay and the transmitters. 
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the relay helps the network to achieve the int-free diversity, e.g., 2 in the (2, 1, 2, 2) MARN (Fig. 5), in 
addition to an improved array gain compared to Scheme 1, e.g., about 10 dB in the (2, 1, 2, 2) MARN 
(Fig. 5) and 2 dB in the (2, 2, 2, 1) MARN (Fig. 6). However, the relay needs to jointly decode both 
user's symbols. Two symbols with 8PSK constellation are jointly recovered in the (2, 1, 2, 2) MARN, 
and four symbols with 8PSK constellation in the (2, 2, 2, 1) MARN. Since the relay does not need to 
decode in Scheme 1, the complexity of Scheme 7 is much higher compared to Scheme 1, especially 
in the (2, 2, 2, 1) MARN. 

VI. Conclusions 

This paper is concerned with multi-user transmission and detection schemes for multi-access relay 
networks, in which multiple sources communicate with one destination by a common relay in two 
hops. For complexity considerations, the nodes in the network have two linear constraints: the relay 
generates its forward signals by linearly transforming its received signals; the destination has linear 
decoding complexity in the number of sources. A new scheme, called IC-Relay-TDMA, is proposed to 
cancel interference at the relay and forward the int-free observations of sources' information in TDMA 
to the destination. IC-Relay-TDMA efficiently allows multi-users to communicate simultaneously in 
the first hop to enhance transmission rate. Through rigorous analysis and simulations, it is shown that 
IC-Relay-TDMA achieves diversity min{J a (i? a — J + l),R a M}. When the number of destination 
antennas is no higher than J a (l — 7^), the maximum int-free diversity R a M is achievable, with a 
higher symbol rate compared to the full-TDMA-DSTC scheme. 

Appendix: Proof of Lemma 1 

First, we show the scenario when the decoding in (13) is symbol-wise, i.e., n = 1, then the scenario 
for multiple symbol joint decoding, i.e., n > 2. 

For n = 1, H n _i = 1, S u = ^ + and = sj 1 ^ from (13). The ML decoding is symbol wise. 
Let a = 7^-, 7 = £~\ and — = Asf\ The pairwise SER of decoding to sj 1 ^ can be 

written as 

P(s? ->• 4 iy ) = ( E Q (JaP Asi 1} \) , (20) 

fki 'dim \ J 

where Q(x) denotes the Gaussian Q function. Note that Q(x) > min < |e~~, ~ > and both 

"I x ^ 1 x ^ (1) flV 

g^e~~ and |e~~ are decreasing functions. Thus, for any e > 0, P(s\ ' — > s\' ) > Pyy < 
e) min ||e~ aPA2e , |^=^= j with A the minimum distance between any two constellation points. Let 
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e = P _1 , we have P(s[ 1} -> ) > P(7 < e)min jie~ QA2 , |^=^ j> where ( is a constant independent 

S v ' 

c 

of e. Thus, diversity can be calculated as 

d= _ Um log P(4 1} -> gg) < Um log P(7 < e) + log C = Um log Pil < e) ^ 

P-s>oo log P e->0+ log e e^0+ log € 

This shows that the diversity is upperbounded by the right-hand-side (RHS) of (17). Next, we show that 
the diversity is also lowerbounded by the RHS of (17). Denote (5 as the maximum distance between 
any two constellation points. For any e > 0, using the Chernoff bound on SER and noticing that 
e - " 2 P7 is a decreasing function with 7, we have 

2 P e 2 /*00 o 2 

PtfU.f^Ee-^^e-^/W^+l e-^ 2 /( 7 )rf7< J P(7< e )+e-^, (21) 

where 7(7) is the probability density function of 7. Let e = P~ n ,0 < n < 1. As P increases, the 
RHS of (21) is dominated by P( 7 < e). The diversity can be lowerbounded as d > \im^ nlo ^' y e <t K 
Since n can be chosen very close to I, the lowerbound and upperbound converge. 

Now, we consider the case of n > 1. Denote Aso 1 " 1 = — . The pairwise error probability of 
decoding a vector of 2 n_1 symbols to in the ML decoding in (13) can be written as 




V 7 /^.a™ V v v ' 



(22) 



Since H^E^H^x -< tr {H^E^H^}!^-!, we have that T < tr {H^E^H^xjH As^ || 2 < 
2 ra_1 tr {H^-lE^H^i}/? 2 with j3 the maximum distance between any two points in all constellations. 
Thus, a lowerbound on (22) can be obtained as P x = Q ( A /a2 n - 1 Ptr {H^E-iH^i}/? 2 ). The 
diversity of p upperbounds that of P (s^ — > j . Recall that p is similar to (20). The diversity of 
the lowerbound can be evaluated by (17) using tr {H*_ 1 S~ 1 H n _i} as the instantaneous normalized 
receive SNR. Thus, the diversity obtained by using 7 = tr {H*_ 1 S~ 1 H„_i} in (17) upperbounds that 
ofP^Us^). 

Next, we find the lowerbound on the diversity of Since the entries of s ( are 

collected from finite constellations satisfying (1), there exists a positive number 9 that lowerbounds all 



hAs K 



(i) 



2 



Z=l:2 Tl - 1 



. Noticing that E u is diagonal from (14). We have T = £ A^ 1 hjAs^ > £ \^9 



2 



-1/32 



tr {H*_ 1 S u 1 H„_ 1 }2^r with A; the l-th diagonal entry of E u . Therefore, the diversity obtained 
by using 7 = tr {H;_ 1 S u 1 H„_i} in (17) lowerbounds that of P (a® ->■ s^. The diversity of 
P (si 1} ->■ si iy j can be calculated by using 7 = tr {H T *_ 1 S u 1 H n _ 1 } in (17). 
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Performance comparison in (J, J a ,R a , M) MARNs. 
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Scheme 


Concurrent Transmission 


Diversity 


Symbol Rate 


Linearity 


IC-Relay-TDMA 


only the user-relay link 


mm{J a (R a -J+l),R a M} 


J + Rn 


Yes 


full-TDMA-DSTC 


none 


R a min{ J a , M} 


Ho 

■m+Ro) 


Yes 


DSTC joint-user ML decoding 


both links 


R a min{ J a ,M} 


2 


No 


DSTC-ICRec 


both links 


mm{Ja,M - 1} 
for J = 2, J a = 1,2, R a = 2 


1 

2 


Yes 
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Fig. 1. Multi-access relay networks. 
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Fig. 3. Block diagram for IC-Relay-TDMA. 




Fig. 4. BER performance IC-Relay-TDMA, under BPSK modulation. 
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Fig. 5. Performance comparison in a (2, 1,2,2) MARN, under 1 bit/user/channel use. 
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Fig. 6. Performance comparison in a (2, 2, 2, 1) MARN, under 1 bit/user/channel use. 
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